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STREAMFLOW  ROUTING  OPTIMIZATION  (OPROUT) 


USERS  MANUAL 


The  Hydrologic  Engineering  Center 


INTRODUCTION 


1.  ORIGIN  OF  PROGRAM 

The  computer  program  "Streamflow  Routing  Optimization"  (OPROUT)  was  developed 
at  the  Hydrologic  Engineering  Center  (HEC)  by  Vernon  R.  Bonner.  The  program 
consists  of  optimization  routines  developed  originally  by  Dr.  Anthony  Slocum  and 
Ramesk  Danekar  of  Anderson-Nichols  &  Company,  Inc.,  as  a  portion  of  HEC  contract, 
"Evaluation  of  Streamflow  Routing  Techniques  with  Special  Emphasis  on  Determining 
Nonlinear  Routing  Criteria,"  October  1975,  (Contract  No.  DACW05-75-C-0027) .  These 
routines  were  linked  with  routing  and  plotting  routines  from  the  computer  program 
HEC-5,  "Simulation  of  Flood  Control  and  Conservation  Systems,"  into  one  program, 
OPROUT,  which  this  manual  describes. 

2.  PURPOSE  OF  PROGRAM 

This  program  can  be  used  to  determine  Modified  Puls  or  Muskingum  routing 
criteria  by  optimizing  parameters  for  a  single  reach  using  observed  upstream  and 
downstream  hydrographs  for  one  to  five  events.  The  main  advantage  to  the  program 
is  that  the  user  can  develop  routing  criteria  from  several  events  by  making  one 
computer  program  execution. 
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DESCRIPTION 


3.  TECHNICAL  APPROACH 

a.  Modified  Puls  Routing  Method.  In  this  routing  method,  outflow  is  a 
uniqiie  function  of  storage.  The  following  relationships  apply: 

At  10  (1) 

where 

I  =  average  inflow  in  the  reach  during  time  interval  At 
0  =  average  outflow  from  the  reach  during  time  interval  At 
AS  =  change  in  storage  during  time  interval  At 


n  - 

v 


(2) 


where  0-, 
respectively 


and  Og  =  reach  outflow  at  staj  t  and  end  of  time  interval, 

s2  ■  $i 

AS  =  J4-L  (3) 


where  S-,  and  S£  =  reach  storage  at  start  and  end  of  time  interval, 
respectively 

substituting  (2)  and  (3)  in  (1)  and  rearranging. 


h.  +  °i  -  ii  Hi  +  r 

At  2  At  "  2  (4) 


Let  SI  =  S  +  Tj- 


SI2 


=  storage  indication  and  At  =  1  then, 
SI-,  +  I  -  Q, 


(5) 


equation  (5)  can  be  solved  using  a  storage  indication-outflow  (SI  vs  0) 
relationship  and  knowing  reach  inflow. 
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"an  be  seen  from  the  described  relationships,  modified  Puls  routing  requires 
a  storage-outflow  relationship  for  the  downstream  end  of  the  routing  reach.  The 
relation  indicates  the  amount  of  flow  occurring  at  the  downstream  end  of  the  reach 
for  a  given  amount  of  water  stored  in  the  reach  between  the  upstream  and  downstream 
ends . 

Starting  with  the  given  information,  an  observed  hydrograph  at  the  upstream 
end  of  the  reach  and  an  observed  hydrograph  at  the  downstream  end  of  the  reach, 
the  problem  is  associated  with  defining  the  proper  storage-outflow  relationship. 

The  proper  relationship  when  used  with  modified  Puls  routing  will  generate  a 
routed  hydrograph  (observed  upstream  hydroyraph  routed  to  the  downstream  end) 
which  is  the  same  as  would  occur  in  the  real  river  reach.  If  there  were  zero 
local  flows,  the  routed  hydrograph  would  be  equal  to  the  observed  downstream 
hydrograph.  In  reality  there  will  be  local  flow  entering  the  reach  during  at 
least  part  of  the  time  period  of  the  observed  hydrographs.  If  this  local  flow 
was  gaged  (all  local  flow  coming  from  a  gaged  tributary),  the  local  flow  could 
be  subtracted  from  the  downstream  observed  hydrograph.  The  result  would  provide 
the  upstream  routed  hydrograph.  Unfortunately,  in  the  majority  of  cases,  the 
local  flow  is  ungaged  and  not  known.  In  this  situation,  the  observed  upstream 
hydrograph  can  be  routed  to  the  downstream  end  and  subtracted  from  the  observed 
downstream  hydrograph  to  generate  the  local  flow  hydrograph.  If  the  routed 
hydrograph  is  reasonable,  no  negative  local  flows  will  be  generated  by  this 
procedure  (i.e.,  for  all  time  periods,  the  routed  hydrograph  is  always  less 
than  or  equal  to  the  downstream  observed  hydrograph).  This  fact  is  used  to 
form  the  objective  function  for  the  optimization  routine.  It  is  stated: 
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Objective  Function  -  Minimize  the  Sum  of  Negative  Local  Flow  in  a  given 

routing  reach  by  proper  defintion  of  the  storage- 
outflow  relationship  and  number  of  routing  subreaches 
using  Modified  Puls  with  upstream  and  downstream 
observed  hydrographs. 

One  form  for  the  storage-outf1*  relationship  is  as  follows: 

S  =  K0M  (6) 

where 

S  =  total  sotrage  in  the  routing  reach 

0  =  reach  outflow 

K,  M  *  are  constants  which  cause  the  function  to  represent  the 

storage-outflow  relationship  of  the  river  reach  under  study 
With  this  function  a  gradient  search  approach  can  be  used  to  achieve  the 
stated  objective  function.  Computationally,  a  beginning  K  and  M  are  selected 
from  the  observed  hydrograph  information  (M  =  1  and  K  equal  travel  time  of 
peak  flow).  The  number  of  subreaches  is  set  to  one.  A  modified  Puls  routing 
is  performed  based  upon  the  selected  S  =  KO^  function.  The  routed  hydrograph 
is  subtracted  from  the  observed  downstream  hydrograph  which  determines  the  sum 
of  the  negative  local  flow.  Either  K  or  M  is  changed  following  a  gradient 
search  procedure  and  a  new  sum  negative  local  flow  ;s  computed  which  can  be 
compared  to  the  last  value.  This  procedure  is  continued  until  the  objective 
function  can  no  longer  be  improved.  The  number  of  routing  subreaches  is 
incremented  by  one  and  the  gradient  search  procedure  is  continued.  The  result 
provides  the  optimum  K,  M  and  number  of  subreaches  for  the  given  routing  reach 
with  the  given  observed  hydrograph.  The  selected  exponent  usually  is  in  the 
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following  range: 


0.6  <  M  <  0.8 

and  the  number  of  subreaches  usually  equals  2. 

A  curve  warping  routine  was  written  to  overcome  the  computational  difficulty 
and  the  inflexibility  of  using  only  S  =  KO^  to  represent  storage-outflow.  This 
routine  is  very  simple  and  surprisingly  effective. 

The  optimization  routine  still  begins  with  a  macro  gradient  search  S  =  KO^ 
to  obtain  a  beginning  storage-outflow  function  which  is  reasonable.  Then  the 
curve  warping  routine  is  called.  The  computed  storage-outflow  relationship  is 
divided  into  17  piece-wise  linear  segments  (this  is  also  used  in  other  HEC 
programs  using  modified  Puls).  This  forms  a  table  of  storage  versus  outflow 
for  equally  spaced  flow  segments  starting  at  zero  and  ending  at  the  peak  flow 
of  the  observed  downstream  hydrograph.  The  number  of  segments  can  be  input 
by  the  user  up  to  a  maximum  of  18  segments. 

The  curve  warping  routine  works  in  the  following  manner.  One  cycle  consists 

of  testing  each  point  in  the  storage-outflow  table  beginning  with  the  lowest 

point  above  zero.  The  flow  is  held  constant  and  the  storage  is  stepped  up, 
down  or  not  changed  depending  if  the  objective  function  is  improved.  Storage  is 
never  stepped  greater  than  the  value  for  the  next  higher  point  in  the  table  or 
less  than  the  next  lower  point  in  the  table  to  keep  the  function  single  valued. 

If  an  improvement  was  made  during  the  last  cycle,  a  new  cycle  is  performed.  This 

is  continued  until  no  more  improvement  can  be  obtained.  The  step  size  is  reduced 

and  the  process  is  repeated.  When  no  more  improvement  can  be  be  made,  the  number 
of  subreaches  is  incremented  by  one  and  the  curve  warp  routine  is  recalled. 

The  final  result  provides  the  optimum  number  of  subreaches  and  the  storage- 
outflow  function.  The  several  additional  degress  of  freedom  in  this  function 
enable  a  solution  which  works  well  for  all  magnitudes  of  flood  events  within 
those  used  to  derive  the  storage-outflow  relationship. 


The  program  also  employs  an  optional  curve  fit  routine  for  smoothing  the 
developed  storage-outflow  function  after  the  curve  warping  routine.  The  program 
also  incorporates  several  other  components  of  the  objective  function  described 
previously  to  improve  the  reasonableness  of  the  local  flow  hydrograph  and  to 
force  the  routed  hydrograph  toward  the  recession  side  of  the  given,  observed 
downstream  hydrograph.  The  total  routing  objective  function  is  defined  as  follows 
MINIMIZE  SUMNL  =  SUM  (SUM!  +  SUM2  +  SUM3)  (7) 

where  the  three  components  are  defined  as, 

1)  For  all  time  periods, 

SUM!  =  SUM  (negative  local  flows)  *  WT1 

2)  For  the  recession  limb  of  the  routed  hydrograph  defined  as  from 
85%  of  the  routed  hydrograph  peak  discharge  to  15%  of  the 
observed  downstream  hydrograph  peak  discharge,  and  where  the 
routed  flow  is  less  than  95%  of  the  observed  downstream  flow 
(indicating  how  early  this  portion  of  the  routed  hydrograph  is), 

SUM2  =  SUM  (routed  hydrograph  minus  observed  downstream 
hydrograph)  *  WT2 

and  3)  For  all  time  periods  following  the  peak  discharge  of  the 
routed  hydrograph  (indicating  how  late  this  portion  of  the 
routed  hydrograph  is), 

SUM3  =  SUM  (negative  local  flows)  *  WT3 

Weighting  of  the  individual  components  of  the  objective  function  (WT1 ,  WT2  and 
WT3)  can  be  input  by  the  user  on  the  »’J  card  as  described  in  Exhibit  2,  Input 
Description.  One  or  both  of  the  second  (SUM2)  and  the  third  (SUM3)  can  be 
eliminated  by  default  values  of  zero  for  WT2  and/or  WT3  thereby  changing  the 
basis  of  the  objective  function.  The  program  output  lists  the  value  (error) 
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of  each  of  the  three  components  and  the  total  (SUM! ,  SUM2,  SUM3  and  SUMNl)  for 
each  iteration  as  NEG  LOCAL,  TOO  EARLY,  TOO  LATE  and  TOTAL,  respectively. 


b.  Muskingum  Routing  Method.  In  this  method,  outflow  from  a  routing  reach 
is  a  linear  function  of  the  sum  of  prism  and  wedge  storage  in  the  reach.  The 
basic  routing  equation  is: 

°n  =  Vn  +  C2In-l  +  Vn-2  ....  (8) 


where: 

0n  =  Ordinate  of  outflow  hydrograph  at  time  n 
In»  etc.  =  Ordinates  of  inflow  hydrograph  at  times  n,  n-1 ,  etc. 

C.| ,  C2,  etc.  =  Routing  coefficients,  as  coefficients  of  inflow 
Equations  used  to  determine  the  coefficients  C-| ,  C^,  etc.,  are  as  follows: 

C1  =  (At  -  2XK)  /  (2  K(l-X)  +  At) 

CC  =  ( ( 2K( 1 -X )  +  At)  -  2 At )/(2K(l-X)  +  At) 

C2  =  C1  •  CC  +  (At  +  2KX)  /  (2K(1-X)  +  At) 

C.  =  C._,  •  CC  for  i  >  2 

where 

At  =  Routing  time  increment 
K  =  Muskingum  routing  parameter  having  units  of  time 
X  =  Muskingum  dimensionless  routing  parameter  between  0  and  .5 
The  program  has  the  capability  of  optimizing  Muskingum  routing  coeffici  •  • 
using  the  same  techniques  and  search  procedures  as  described  for  the  modified  Puls 
method.  The  M  in  equation  (6)  is  equal  to  1  for  the  linear  Muskingum  method  and 
the  same  objective  function  is  used  to  optimize  Muskingum  routing  constants  and 
K  and  X. 


(9) 

(10) 

(ID 

(12) 
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From  the  above  relationships  it  can  be  seen  that  the  following  relationship 
between  K  for  each  subreach  and  At  must  be  true  to  avoid  negative  coefficients. 


2( 1 -X)  "  At  "  2X 

If  a  Muskingum  routing  optimization  produces  negative  coefficients,  the  user 
should  increase  the  number  of  subreaches,  thereby  reducing  K,  so  that  the  above 
limits  are  met.  The  method  for  defining  the  number  of  subreaches  (RT.3,  3rd 
field  of  the  RT  Card)  is  described  under  paragraph  4,  Program  Capabilities  and 
illustrated  in  Exhibit  1. 
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4.  PROGRAM  CAPABILITIES 

The  program  will  provide  solutions  to  the  following  problems: 

a.  Determine  Muskingum  K  and  X  and  the  number  of  routing  subreaches  for  a 
single  reach  with  up  to  five  setr  of  given  upstream  and  downstream  observed 
hydrographs.  Each  hydrograph  set  is  given  equal  weight  in  the  optimization 
objective  function. 

b.  Determine  Modified  Puls  storage-outflow  relationship  and  number  of 
routing  subreaches  for  a  given  reach  with  up  to  five  multiple  sets  of  given 
upstream  and  downstream  observed  hydrographs. 

c.  Determine  routing  criteria  for  a  given  number  of  subreaches  in  either 
cases  a  or  b  above. 

d.  Determine  the  storage-outflow  curve,  as  in  case  b  above,  with  a  given 
coefficient  X  (Working  R&D  routing)  or  an  additional  lag  of  the  routed  hydrograph. 

Under  any  of  the  above  ortions,  a  complete  trace  feature  is  available  to 
monitor  the  progress  of  the  optimization  computations.  This  feature  is  written 
in  the  same  manner  as  exists  in  HEC-5;  therefore,  a  source  listing  is  required 
to  interpret  results. 

The  program  can  be  used  to  solve  the  more  difficult  problem  when  two  or  more 
upstream  gages  and  routing  reaches  flow  to  a  common  downstream  gage.  This 
currently  cannot  be  solved  automatically  by  the  program;  however,  the  user  can 
develop  routing  criteria  for  this  situation  by  multiple  executions  of  the  current 
program.  In  most  cases,  one  upstream  gage  will  have  the  dominant  flow.  The 
optimization  routine  can  be  operated  using  the  dominant  upstream  gage  as 
the  observed  upstream  hydrograph  and  the  given  observed  downstream  hydrograph. 

The  resulting  local  flow  hydrograph  which  is  computed  by  the  Optimization  Routine 
can  be  used  as  the  observed  downstream  hydrograph  in  conjunction  with  the  next 
largest  upstream  gage  flow  which  becomes  the  observed  upstream  hydrograph. 
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This  process  is  repeated  until  routing  criteria  is  defined  for  all  routing  reaches. 
In  most  cases,  some  user  smoothing  of  the  local  flow  hydrographs  will  be  required 
(while  maintaining  consistent  volume)  and  a  second  iteration  will  be  required 
beginning  with  the  dominant  upstream  gage.  The  problem  is  difficult  to  solve 
due  to  the  many  additional  degress  of  freedom.  For  each  additional  routing 
reach,  two  basic  unknowns  are  added,  routing  criteria  and  local  flow;  whereas, 
only  one  known  is  added,  upstream  observed  hydrograph. 

During  the  performance  of  the  optimization  routine,  a  printer  plot  can  be 
requested  (5th  field  of  J1  card)  which  will  plot  the  upstream  observed  hydrograph 
(values  entered  on  IN  cards),  downstream  observed  hydrograph  (values  entered  on 
IN  cards),  routed  hydrograph  (upstream  hydrograph  routed  to  downstream  location 
using  routing  criteria  derived  by  optimization  routine)  and  the  computed  local 
flow  hydrograph  (difference  between  routed  hydrograph  and  downstream  observed 
hydrograph).  This  provides  a  visual  check  on  the  results  developed  by  the 
optimization  subroutine. 

In  addition,  upon  completion,  the  optimization  routine  prints  the  values  for 
the  adopted  storage-outflow  relationship  (table  of  storage  versus  flow)  and  the 
optimum  number  of  subreaches.  The  travel  time  indicator  (inverse  slope  of 
storage-outflow  relationship)  is  also  printed  for  each  incremental  linear  segment 
of  the  storage-outflow  table  and  for  the  given  point  to  the  origin.  The  upstream 
hydrograph  is  printed  as  read  in  addition  to  the  routed  hydrograph.  Incremental 
local  flows  are  printed  twice  for  each  control  point.  The  first  is  the  computed 
values  without  adjustment  and  the  second  is  the  adjusted  values  with  all  computed 
negative  values  set  to  zero.  The  negative  volume  is  proportioned  to  the 
remaining  positive  values.  Program  capabilities  are  illustrated  in  Exhibit  1  - 
Example  Input  and  Output.  A  detailed  description  of  the  program  input  is 
presented  in  Exhibit  2. 
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EXHIBIT  1 


EXAMPLE  INPUT  AND  OUTPUT 


EXHIBIT  1 


EXAMPLE  INPUT  AND  OUTPUT 
Streamflow  Routing  Optimization  (OPROUT) 

The  input  and  output  for  three  examples  are  provided  to  illustrate  the  use 
of  selected  program  capabilities  and  options  and  to  assist  in  verifying  the 
correct  execution  of  the  program.  A  brief  description  of  each  example  is 
provided  below.  Printer  plots  are  requested  for  each  of  the  examples. 

a.  Example  1  -  To  determine  modified  Puls  storage-outflow  relationship 
and  the  number  of  routing  subreaches  for  a  single  reach  and  from  a  single 
flood  event. 

b.  Example  2  -  To  determine  modified  Puls  storage-outflow  relationship 
and  the  number  of  routing  subreaches  using  three  flood  events.  (J1 .7) 

c.  Example  3  -  To  determine  Muskingum  K  and  X  for  one  flood  event  and  a 
user  specified  number  of  routing  subreaches.  (RT.3) 


Exhibit  1 
Example  1 
1  of  10 


EXAMPLE  1  INPUT 


TlEXAMPLE  1  MODIFIED  PULS  ROUTING  OPTIMIZATION 
T2  SCIOTO  RIVER  CHILLICOTHE  TO  HIGBY  OHIO 
T3  SINGLE  EVENT  -  JANUARY  1959 


J1  100 

3 

0 

0 

1 

0 

0 

0 

1 

1 

ID  CHILLICOTHE  TO 

HIGBY  - 

JAN  1959 

RT  2 

3 

.5 

0 

0 

0 

6 

3 

3 

0 

IN  2 

JAN  59 

QO  3160 

3093 

3025 

2940 

3020 

3207 

3800 

4660 

5650 

7203 

QO  8395 

9653 

11400 

13300 

15200 

17075 

18600 

20150 

21700 

23575 

QO  25800 

31250 

42800 

64650 

99900 

126000 

140000 

143000 

138000 

130500 

Q0120500 

110500 

103000 

94300 

85500 

76950 

69000 

63100 

57200 

52350 

QO  47500 

43950 

40400 

37300 

34200 

31750 

29300 

27550 

25800 

24450 

QO  23100 

22000 

20900 

19950 

19000 

18300 

17600 

17150 

16700 

16500 

QO  16300 

16225 

16150 

16075 

16000 

15900 

15800 

15650 

15500 

15300 

QO  15100 

14800 

14500 

14125 

13750 

13375 

13000 

12725 

12450 

12175 

QO  11900 

11600 

11300 

11000 

10750 

10500 

10250 

10000 

9750 

9500 

QO  9250 

9000 

8750 

8500 

8250 

8000 

7750 

7500 

7250 

7000 

IN  3 

JAN  59 

QO  4580 

4449 

4317 

4277 

4420 

5850 

7635 

9595 

11800 

14500 

QO  17750 

21450 

25200 

28125 

31550 

35000 

37900 

39580 

41450 

42975 

QO  44700 

47050 

49400 

53300 

59400 

74850 

97350 

121000 

143000 

156500 

Q0160000 

157000 

149000 

138000 

127000 

115000 

103000 

93200 

83400 

76850 

QO  70300 

64950 

59600 

55550 

51500 

47750 

44000 

40750 

37500 

34950 

QO  32400 

29850 

27300 

26000 

24700 

23400 

22100 

21300 

20500 

20000 

QO  19500 

19200 

19025 

18912 

18800 

18300 

18300 

18300 

18300 

18300 

QO  18300 

18300 

18300 

15700 

15700 

15700 

15700 

15700 

15700 

15700 

QO  15700 

13100 

13100 

13100 

13100 

13100 

13100 

13100 

13100 

11200 

QO  11200 

11200 

11200 

11200 

11200 

11200 

11200 

8600 

8600 

8600 

ROUTING  OPTIMIZATION  INPUT  DATA 
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EXAMPLE  2  INPUT 


TlEX AMPLE  2  MODIFIED  PULS  ROUTING  OPTIMIZATION 


T2  SCIOTO  RIVER 

CHILLICOTHE  TO 

HIGBY  OHIO 

T3  THREE 

EVENTS 

-  MAR  1964,  FEB 

1959,  JAN  1959 

J1  100 

3 

0 

0 

1 

0 

3 

0 

1 

1 

ID  CHILLICOTHE  TO 

HIGBY  - 

3  EVENTS 

RT  2 

3 

.5 

0 

0 

0 

6 

3 

3 

0 

IN  2 

MAR64 

QO  333 

330 

400 

500 

500 

550 

750 

1250 

2250 

4000 

QO  6500 

9000 

11500 

13000 

14000 

14500 

15000 

15500 

15750 

16000 

QO  16500 

16750 

17000 

17000 

17000 

17000 

17000 

16800 

17000 

17000 

on  17000 

17000 

16700 

16250 

15500 

15000 

14000 

13250 

12000 

11200 

Qo  10500 

10200 

10500 

11700 

14200 

17000 

19500 

21500 

23700 

26500 

QO  30000 

34000 

40000 

45500 

50000 

53000 

55200 

58300 

61500 

65000 

QO  68500 

70750 

71600 

71500 

70600 

69000 

66000 

62500 

59000 

55500 

QO  51750 

48000 

45200 

42000 

39500 

37000 

33750 

30750 

28500 

26000 

QO  22600 

22100 

23000 

22500 

21500 

20750 

20000 

19500 

19000 

18750 

QO  18500 

18000 

17500 

17000 

16500 

16250 

16000 

16000 

16000 

16000 

IN  3 

MAR64 

QO  628 

633 

637 

993 

1350 

3787 

6225 

8663 

11100 

15350 

QO  19600 

22450 

25300 

27600 

29900 

31150 

32400 

31700 

31000 

30300 

QO  29600 

28900 

28200 

27500 

26800 

26300 

25800 

25300 

24800 

24275 

QO  23750 

23225 

22700 

21850 

21000 

20150 

19300 

18375 

17450 

16252 

QO  15600 

15650 

15700 

22000 

28300 

35250 

42200 

47450 

52700 

59100 

QO  65500 

72350 

79200 

91600 

104000 

115000 

126000 

125000 

125000 

124500 

Q0124000 

121000 

118000 

115000 

112000 

109500 

107000 

102850 

98700 

93375 

QO  88050 

82725 

77400 

72575 

67750 

62925 

58100 

54575 

51050 

47525 

QO  44000 

41775 

39550 

37325 

35100 

34200 

33300 

32400 

31500 

31000 

QO  30500 

30000 

29500 

28625 

27750 

26875 

26000 

25425 

24850 

24275 

IN  2 

FEB59 

QO  3310 

3310 

3310 

3310 

3310 

3310 

3310 

3310 

3310 

3360 

QO  3480 

4580 

5680 

7670 

9660 

11630 

13600 

14600 

15600 

16750 

QO  17900 

19675 

21450 

23225 

25000 

28400 

31800 

35200 

38600 

40700 

QO  42800 

43100 

42600 

40575 

38550 

36525 

34500 

32300 

30100 

28350 

QO  26600 

25350 

24100 

23200 

223r  ? 

22125 

21950 

21775 

21600 

21500 

QO  21400 

21300 

21200 

21175 

21150 

21125 

21100 

21300 

21500 

21700 

QO  21900 

22000 

22100 

22000 

21900 

21650 

21400 

20850 

20300 

19600 

QO  18900 

18200 

17500 

17013 

16525 

16037 

15550 

15063 

14575 

14087 

QO  13600 

13275 

12950 

12625 

12300 

11850 

11400 

10620 

9840 

8860 

QO  7880 

7215 

6550 

6283 

6015 

5747 

5500 

5300 

5100 

4900 

IN  3 

FEB59 

QO  5210 

5210 

5210 

5210 

5210 

5210 

5210 

5210 

5210 

5480 

QO  5890 

6300 

8050 

13175 

17300 

20700 

24100 

25450 

26800 

26600 

QO  26400 

26200 

26000 

27000 

28000 

29700 

31400 

33850 

36300 

38150 

QO  40000 

41850 

43700 

43850 

44000 

43400 

42800 

41500 

40200 

38450 

QO  36700 

34900 

33100 

31650 

30200 

29550 

28900 

30100 

31300 

32200 

QO  33100 

32250 

31400 

30550 

29700 

28850 

28000 

27650 

27300 

27250 

QO  27200 

27150 

27100 

26950 

26800 

25000 

25000 

25000 

25000 

25000 

QO  25000 

25000 

25000 

23937 

22875 

21813 

20750 

19687 

18625 

17563 

QO  16500 

15875 

15250 

14625 

14000 

13375 

12750 

12125 

11500 

11063 

QO  10625 

10187 

9750 

9313 

8875 

8437 

8000 

7600 

7200 

6800 
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IN  t  JAN59 


QO  3160 

3093 

3025 

2940 

3020 

3207 

3800 

4660 

5650 

7203 

QO  8395 

9653 

11400 

13300 

15200 

17075 

18600 

20150 

21700 

23575 

QO  25800 

31250 

42800 

64650 

99900 

126000 

140000 

143000 

138000 

130500 

Q0120500 

110500 

103000 

94300 

85500 

76950 

69000 

63100 

57200 

52350 

QO  47500 

43950 

40400 

37300 

34200 

31750 

29300 

27550 

25800 

24450 

QO  23100 

22000 

20900 

19950 

19000 

18300 

17600 

17150 

16700 

16500 

QO  16300 

16225 

16150 

16075 

16000 

15900 

15800 

15650 

15500 

15300 

QO  15100 

14800 

14500 

14125 

13750 

13375 

13000 

12725 

.  12450 

12175 

QO  11900 

11600 

11300 

11000 

10750 

10500 

10250 

10000 

9750 

9500 

QO  9250 

9000 

8750 

8500 

8250 

8000 

7750 

7500 

7250 

7000 

IN  3 

JAN59 

QO  4580 

4449 

4317 

4277 

4420 

5850 

7635 

9595 

11800 

14500 

QO  17750 

21450 

25200 

28125 

31550 

35000 

37900 

39580 

41450 

42975 

QO  44700 

47050 

49400 

53300 

59400 

74850 

97350 

121000 

143000 

156500 

QO160000 

157000 

149000 

138000 

127000 

115000 

103000 

93200 

83400 

76850 

QO  70300 

64950 

59600 

55550 

51500 

47750 

44000 

40750 

37500 

34950 

QO  32400 

29850 

27300 

26000 

24700 

23400 

22100 

21300 

20500 

20000 

QO  19500 

19200 

19025 

18912 

18800 

18300 

18300 

18300 

18300 

18300 

QO  18300 

18300 

18300 

15700 

15700 

15700 

15700 

15700 

15700 

15700 

QO  15700 

13100 

13100 

13100 

13100 

13100 

13100 

13100 

13100 

11200 

QO  11200 

11200 

11200 

11200 

11200 

11200 

11200 

8600 

8600 

8600 
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EXAMPLE  2  OUTPUT 
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oooooooooo 


oooooooooo 

ooooooooino 

ooocMtnomor-o 

■^•lor-iHvDininvDoovr) 


oomcNJootntnotn 

mot-moor-cMot"* 

nnNNHinninoN 
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ooooooor*oo 
iDinotnooocoioo 
coc-r'-coinr-oocoo 
nioocoHHcn^co^ 
H  N  N  N  H  H 


oooooooooo 

oooooooooo 


oooooooooo 

minoooooooo 

cNr-oor-inoinoo 

cMinr-cMmrHo\comio 


oooooooooo 

ooomoooinom 

rHoco^r-or-omoo 

HH^hNinCOHHV 

HnoiriirtNininnN 


oooooooinoo 
r-toooooor-^ro 
nvo  Hiow  nw  oo  h 
ninoooHHo^min 

H  fO  MN  N  NH 


oooooooooo 

oooooooooo 


oooooooooo 

inootnooomoo 

csincoNinntnMno 

Hin^mHCONOffli£l 

HHHNiniOnHH 


coooinoooinom 

icoor-momr-ocM 

«>r*nn'd'ooointr^ 

00HtnO3Mn(S<fNin 

nNHitNOinfON 


ooootnoorooo 

HOOOhOWWNO 

fOvDNrnhncoovon 

n^mNHHoinoin 

HfOfONfNNrlH 


oooooooooo 

inooooooinoo 

t'-oooincMor-oo 

inr',«r<j\m\ocooi£> 

HrlHHUMOfONH 


tnooooooooo 

CMOOOOOOOOO 
N-<f  OOniMOOHMO 

vocsincncNvor'OOfOU) 

MNH^NOUIMOJ 


OOOOOOOOOO 

.-loootnootnoo 

mmiHTT^Wr^unr-iin 

rlnntSMOJHH 


oooooooooo 

inoooooooinin 

ininoooooor»CM 

^r^inr^romr^ovo 

HHHHt/MOfONH 


r^ooooooinom 

oomominooNoh 

t^rlMHOJOIflO\NCO 

corHiDotntncftCMrrio 

noiMroHovoncM 


ooomininor-'or-« 

rinoNNNinninir 

CO«)^finrHiHU>0000 

nHCO^NHH^Hin 

HNnOlNNHrl 


o  o  o  o  o  o  o 

o  o  o  o  o  o  o 

tO  O  O  10  CM  O  VO 

I'Mn^oo 
rl  H  H  H  in  Is 


OOOOOOOOOO 

inoooooomoin 

fOOOOOPOOOh'rHh- 

HovnrHoo^rtNir-inr^ 
N(.NNOH«5fON 
rH  rH 


ooooooomom 
HwomoinonoH 
nvDOinnHounno 
fnaunOONHHVCNW 
CM  ro  CM  CM  CM  rH  fH 


o  o  o  o  o  o  o  * 

o  o  o  in  o  o  o  < 

mooNMnini 
ro  r-  id  H  I/)  H  I 

rl  rl  rl  H  ’C  r*  1 


roooooooininLn 

o^ooLnooor-cNjcM 

cnioincooiooincoio 

r-r-rHCMrHincMr-'Co 

oi(N>cNfNcn^Hr>roc\j 


ooinmoinowinn 
.Hr-cMr-or-orHcMoo 
fOWNinciHOOWN 
cor-roororHCMr-cMvo 
N  tTN  CJ  NH  H 


o  o  o  O  O  O  o 

o  o  o  o  o  o  o 

t**  in  o  ko  cm  o  in 

ooo*Hir>ror- 

H  rl  V  Is  'T  N  H 


r-000000000 

rooooooooino 

ococMr-r-cMo^rinm 

incoNinwcohO'iOv 

CMCMCMrHt'-rHr-cocM 


oooooooooo 

rHcoinooooomin 

rOVD^kPHCNHlDO^lD 

roinrHCM^iHCMr-cMio 

N  ^  N  N  N  rt  rl 


o  o  ©  o  o  o  o 

o  o  o  in  o  o  o 

o  cm  o  r-  o  h  © 

r-  o  -m*  o  co  cm  co 

HHMhM'NH 


moomoooininoi 

IminocMinmocMr-oi 

lo^mcMvocoot^r^o 

CMOOCOinCMrHCMfHO 

CMCMCMiHr-CMCO^rCO 


looinoooooinin 

iHcor-oinooor-H 

n^vDHnnocNMN 

co«^<r>nm<-tfMcopor' 

rH  CM  CM  CM  rH  H 


•  CM  o  o  O 

CM  O  O  O 


OOOOOOOCOOOOOOOOOOOCMOOOOOOOOOO 
OOOOOOO  OOOOOOOOOO  OOOOOOOOOO 


COOOOOOOOOO 

coooooooinoo 

ooininotnotnr-ioin 

vOOC-OOCOrlCNCO 

HHHPMOU1NH 


coooooooooo 
cMoomoooinoo 
i£>inior*i£>tnooom 
(nmnuunMoovo 
rH  CM  CM  rH  H>  CM  CO  V  n 


oooooooooo 

(HCOOOOOOOOCO 

n’Tcr,  cc'o^cnoooo 
rocor-CMiDrHrHcocnr- 
rH  CM  CM  CM  H  rH 
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ooooooooo 

ooooooooo 


oooooococno 

ointnoinouiuio 
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vooococjr-inr-Hvo 
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cm  iniDrocinnHino 
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OOOOOOOOOO 
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rHrHcocnp-oooincoco 
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oooooooo 


oooooor>in 

ininooinocoN 

TTOOlflHVOOlDH 

innHOMno\N 

Nn^friNNHH 


o  o  o  o  o  o 

o  o  o  o  o  o 


o  o  o  o  o  o 

o  unnoo  in 

ID  IDHOH  in 

r"  -<r  o  mm  p 

CM  VO  CM 


o  o  o  o  o  CO 

O  o  O  O  O  O  O 


o  o  in  o  o  in  o 

in  in  cm  o  o  m  co 

h  vo  p-  o  in  in  in 

o*  in  cm  o  0*  cn  cr> 

rH  i— I  rH  rH  CO 


OOOOOOOO 

OOOOOOOO 


OOOOOOOO 

ootnooooo 
OCMO*fOCnO*rHlO 
rinoHoomnco 
(M  <71  N  H  i-l  H 


ooooooooo 

ooooooooo 


ooooooooo 

ooooooinino 

H^rooinooppo 

-^•rHCMOOCOinOCMOO 

(NCn^TCMCMCMCMH 


ooooooooo 

ooooooooo 


ooooooroinp' 

oooininoHhn 

pp^inooooonvr 

ouinoicoinHfnco 

CINM-CJNOJOiH 


ooooooooo 

ooooooooo 


oooooonom 

oooooor-or^ 

moocMr-cooooco 

ocDTroovvocM^co 

HNvTMNCMNH 


oooooooooo 

oooooooooo 


oooooooooo 

ooooooooinin 

COVDOOmiOCOOCMt^ 

ro  oo  o  OHJi  Mn  n  o  h 

H^yorMHHHrH 


oooooooooo 

oooooooooo 


p'tnooooomoo 

or-oininoor-oo 

(Moooipnmnmo 

n  mo  w  h  to  in  n  o  co 

HNpfOHHHH 


oooooooooo 

oooooooooo 


oooooooooo 

CMOOOOOOininin 

ONOUONOOPPN 

n  in  cwn 'i*  ciio  n  o  to 

HClflOfOrlrlHH 


oooooooooo 

oooooooooo 


inooooooooo 

momooooooo 

U)0\fOOOrlDPHCM 

pppn^oJcoinmH 

fnffiO'fNHrlHH 


oooooooooo 

oooooooooo 


oooooooooo 
inoinoinooooo 
OOOlDOhvf  P1PHN 
inmvrinhmcoinrnH 

H 

oooooooooo 

oooooooooo 


oooooooooo 

cMtnoooooooo 

irinvroinpcoPHN 

vrHffipHtfcotnpiH 

fOinCJlDOJHHHH 


oooooooo 

oooooooo 


inooooor*m 

p-oininminmcM 

Horaicincnouo 

cnr-mHoioco^r 

HCJvrnnNNH 


o  o  o  o  o  o 

o  o  o  o  o  o 


CO  o  o  o  o  o 

iH  Tf  o  m  o  o 

co  c\  co  vo  co  co 

<J\  CM  CO  *T  r- 

rMOOl  CO 


o  o  o  o  o  o  o 

o  o  o  o  o  o  o 


o  in  in  o  o  r-  in 

in  r-  cm  o  o  r-  cm 

OlO  H  Oin  CM  rH 

OMO  M1  H  CO  KT  00 

»— I  rH  rH  rH  CM 


oooooooo 

oooooooo 


OOOOCMOOO 

oonorHooo 

moinoa'PHCS 

coooinvnooincoH 

innunNHHHH 


oooooooooo 

oooooooooo 


ooooooooo 

moooooointn 

OOPH'frlONfs 

covomcnHhininov 

CM  Tf  CO  m  CM  CM  H 


inooooooooo 

iMoooootnoom 

o^coo^mHinrop 

(OHNcnOOlO^H© 

HM'OM'NHHH 


r-ooooomooo 

rHOOOOOCMOOO 

cncM^owmonHN 

■^rincncnor-ovcocnrH 

N^’tinCMHHHH 

rH 


ooooooooocnoooooooooocnoooooooooo 

oooooooooinoooooooooomoooooooooo 

oooooooinr*'  §cncnooooinooo§avooooooooo 
ooinoinmoMDbCHninoinocMooo  ►o^minoinmoooo 

fONffl^NrlOCOH  OVCCMinOVOCMCOVOO  vrvfoociooNcnHN 

coioH^OpiAino  naiHoonNiO'fHci  «h 

CMrrcnmCMCMrHH  fO  H  *t  CM  H  H  H  N^inVf)(MHHHH 


OOOOOOOOOCMOOOr'  OOOOOOCOOOOOOOOOOO 
OOOOOOOOO  OOOOOOOOOO  OOOOOOOOOO 


ooooooooin 

OVOOOOOOOCM 

co^opHNOinic 

in  io  o  in  cn  Mn  vo  o 

CMTTCncOCMCMrHrH 


omoooooooo 
vdctioooooooio 
HcooDininHnHcncM 
ncoinopruDinHtn 
CM  <M  CM  rH  rH  rH 


OOOOOOOOOO 

comoooooooo 

inpponM’incopN 

'Th'roo(M<n®mH 

HVVOPCOHHHH 
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EXAMPLE  3  INPUT 


TlEXAMPLE 

3  MUSKINGUM 

ROUTING 

OPTIMIZATION  USER  SPECIFIED  NO. 

OF  SUBREACHES 

T2  SCIOTO  RIVER 

CHILLICOTHE  TO 

HIG.  Y  OHIO 

T3  SINGLE  EVENT 

-  JANUARY  1959 

J1  100 

3 

0 

0 

1 

0 

0 

0 

1 

1 

ID  CHILLICOTHE  TO 

HIGBY  - 

JAN  1959 

RT  2 

3 

3.6 

0 

0 

0 

6 

3 

3 

0 

IN  2 

JAN59 

QO  3160 

3093 

3025 

2940 

3020 

3207 

3800 

4660 

5650 

7203 

QO  8395 

9653 

11400 

13300 

15200 

17075 

18600 

20150 

21700 

23575 

QO  25800 

31250 

42800 

64650 

99900 

126000 

140000 

143000 

138000 

130500 

Q0120500 

110500 

103000 

94300 

85500 

76950 

69000 

63100 

57200 

52350 

QO  47500 

43950 

40400 

37300 

34200 

31750 

29300 

27550 

25800 

24450 

QO  23100 

22000 

20900 

19950 

19000 

18300 

17600 

17150 

16700 

16500 

QO  16300 

16225 

16150 

16075 

16000 

15900 

15800 

15650 

15500 

15300 

QO  15100 

14800 

14500 

14125 

13750 

13375 

13000 

12725 

12450 

12175 

QO  11900 

11600 

11300 

11000 

10750 

10500 

10250 

10000 

9750 

9500 

QO  9250 

9000 

8750 

8500 

8250 

8000 

7750 

7500 

7250 

7000 

IN  3 

JAN59 

QO  4580 

4449 

4317 

4277 

4420 

5850 

7635 

9595 

11800 

14500 

QO  17750 

21450 

25200 

28125 

31550 

35000 

37900 

39580 

41450 

42975 

QO  44700 

47050 

49400 

53300 

59400 

74850 

97350 

121000 

143000 

156500 

Q0160000 

157000 

149000 

138000 

127000 

115000 

103000 

93200 

83400 

76850 

QO  70300 

64950 

59600 

55550 

51500 

47750 

44000 

40750 

37500 

34950 

QO  32400 

29850 

27300 

26000 

24700 

23400 

22100 

21300 

20500 

20000 

QO  19500 

19200 

19025 

18912 

18800 

18300 

18300 

18300 

18300 

18300 

QO  18300 

18300 

18300 

15700 

15700 

15700 

15700 

15700 

15700 

15700 

QO  15700 

13100 

13100 

13100 

13100 

13100 

13100 

13100 

13100 

11200 

QO  11200 

11200 

11200 

11200 

11200 

11200 

11200 

8600 

8600 

8600 
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EXHIBIT  2 


INPUT  DESCRIPTION 

Streamflow  Routing  Optimization  (OPROUT) 


To  determine  the  optimization  routing  criteria  using  modified  Puls  or 
Muskingum  routing  methods  for  a  single  reach,  at  least  one  set  of  observed 
hydrographs  are  required  for  the  upstream  and  downstream  ends  of  the  reach. 

A  maximum  of  five  sets  of  observed  hydrographs  can  be  used  to  develop  the 
routing  criteria.  The  observed  upstream  hydrograph  must  be  input  ahead  of 
the  downstream  hydrograph  for  each  event.  There  is  a  date  variable  on  the 
input  cards  so  that  the  starting  date  of  each  flood  can  be  input  with  the 
flow  data. 

If  several  floods  are  used  to  develop  the  routing  criteria,  all  the  flow 
data  must  have  the  same  number  of  periods.  This  means  shorter  events  must 
be  extended  with  base  flow  values  up  to  the  number  of  periods  of  the  longest 
event.  The  added  data  should  be  included  at  the  end  of  the  actual  flood  and 
the  upstream  values  generally  should  be  less  than  the  downstream  values. 

The  following  INPUT  DESCRIPTION  provides  a  detailed  description  of  data 
input.  Input  format  is  the  standard  HEC  ten  fields  of  8  columns  each  with 
the  first  two  card  columns  reserved  for  the  card  identification.  Card  ID 
is  not  read  by  the  program.  Three  title  cards  are  used  to  provide  labeling 
information.  The  job  card  defines  job  options  using  integers  (no  decimal 
points).  The  identification  card  provides  a  label  for  plots,  and  the  routing 
card  provides  the  routing  optimization  information.  The  remaining  cards  contai 
the  flow  data  (upstream  then  downstream  hydrograph). 
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INPUT  DESCRIPTION 


T1 

T2 

T3 

J1 


Streamflow  Routing  Optimization  Program  (OPROUT) 

Tl,  T2,  T3  Cards  -  Title  Cards 

Job  title  cards;  three  cards  required.  Both  alphabetic  and  numberic 
information  may  be  placed  on  these  cards.  Information  on  these  cards  will  be 
printed  out  as  job  title  on  the  first  page  of  output. 


J1  Card  -  Job  Card  - 

Field  Variable 

1  NPER 

2  IPER 

3  INPUT 


all  integers 

Value  Description 

+  Number  of  values  on  each  set  of  flow  data 

cards  (IN).  Each  set  of  flow  data  must  have 
NPER  Values  (300  maximum). 

+  Time  interval  (At)  in  hours  between  flow  values. 

0  First  input  card  for  each  flow  set  will  have 

control  point  number  in  field  1  and  date  in 
field  2  and  no  other  data.  Subsequent  cards 
contain  10  fields  of  flow  data. 

5  Flow  data  will  be  in  default  HEC-5  format. 

First  two  fields  of  first  card  are  for  control 
point  and  date  with  flow  data  on  the  remaining 
8  fields.  Subsequent  cards  contain  10  fields 
of  flow  data. 
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J1 


J1  Card  - 
Field 
4 


5 


6 

7 


8 


9 


Continued 


Variable  Value 
IPRNT  0 


1 

15 

20 

I  PLOT  0 
1 


IPUNCH  0 

1 

NFLOOD  0 

+ 

NPTSQ  0 


ICURV  0 


Description 

Output  will  be  limited  to  Input  Data,  Optimization 
Errors,  and  optimization  results. 

Same  as  above  plus  optimization  trace. 

Same  as  above  plus  routing  trace. 

Same  as  above  plus  detailed  routing  trace. 

No  hydrograph  plot. 

Plot  hydrograph  for  Inflow  (Routed  Downstream), 
Observed  flow  downstream,  Computed  Local  Flow, 
and  Inflow  Upstream. 

No  punch  cards. 

Punch  computed  local  flow  data. 

One  flood  event  will  be  used  to  optimize  routing 
criteria. 

Number  of  flood  events  (5  maximum)  requires 
upstream  and  downstream  flow  data  for  each  event. 
Number  of  storage-outflow  points  to  be  derived 
from  optimization  routing  is  equal  to  9. 

Number  of  storage-outflow  points  to  be  derived 
(18  maximum).  The  more  points  used,  the  longer 
the  optimization  time. 

No  adjustments  will  be  made  to  the  storage-outflow 
curve  (UNSMOOTHED  CURVE). 
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J1  Card  -  Continued 


J1 

ID 

RT 


Field  Variable  Value 

+ 


10  I FLOW  0 

1 


Description 

A  fourth  order  polynomial  curve  should  be 
fitted  to  the  unsmoothed  curve.  The  smooth 
curve  will  then  be  used  for  computing  local 
fl ows. 

Input  flow  data  is  average  for  the  period. 
Input  flow  data  is  end  of  period  data. 

Flow  data  will  be  averaged  before  routing. 


ID  Card  -  Identification  Card  for  Reach 

Field  Variable  Value  Description 

1-  5  CPT(K)  +  Title  (alphanumeric)  of  routing  reach  in 

columns  3-40.  Title  will  be  printed  on 
hydrograph  plots. 

6-10  Not  used 


RT  Card  -  Routing  Card 

Field  Variable  Value 

1  RTFR  + 


2  RTTO  + 


3  RTMD  + 


Description 

Control  point  number  of  upstream  end  of  routing 
reach. 

Control  point  number  of  downstream  end  of 
routing  reach. 

Number  of  subreaches  to  the  left  of  the 
decimal  and  routing  method  to  the  right  of 
the  decimal  (.5  for  modified  Puls  and  .6  for 
Muskingum).  If  subreaches  are  not  given, 
the  number  will  be  optimized. 


4 


* 
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RT 


RT  Card  - 

Continued 

Field 

Variable 

Value 

Description 

4 

RTCOF 

0 

+ 

Muskingum  routing  coefficient  "X"  will  be 

optimized  if  Muskingum  optimization  is  requested 

(RT.3  =  .6). 

Coefficient  X  will  not  be  optimized.  The  given 

value  will  be  used  for  Muskingum  or  modified  Puls 

5 

XMUSK 

0 

Travel  time  (Muskingum  K)  in  hours  will  be 

optimized  if  Muskingum  optimization  is  requested. 

6 

LAG 

0 

+ 

No  lag  in  addition  to  routing. 

In  addition  to  routing,  lag  outflow  by  the 

number  of  periods  shown. 

7 

WT1 

0 

+ 

Weighting  of  negative  local  flow  equals  one  (1.0) 

Weighting  factor  for  negative  local  flow. 

Use  any  number. 

8 

WT2 

0 

+ 

Weighting  factor  (penalty)  for  recession  leg 

being  too  early  equals  zero  (0). 

Weighting  factor  for  recession  leg  being  too 

early.  (Suggested  value  near  1.0.) 

9 

WT3 

0 

+ 

Weighting  factor  of  error  on  recession  leg  of 

hydrograph  being  too  late  equals  zero  (0). 

Weighting  factor  for  recession  leg  being  too 

late.  Use  in  conjunction  with  WT2.  (Suggested 

value  near  4.0.) 

10 

METRIC 

0 

+ 

English  units. 

Metric  units. 

Exhibit  2 
5  of  6 


IN 


IN  Cards  -  Flow  data  for  upstream  and  then  downstream  station.* 


Field 

Variable 

Value 

Description 

1 

MM 

+ 

Control  point  number  for  input  hydrograph. 

First  hydrograph  for  the  upstream  location, 

and  then  the  downstream  location. 

2 

DATE 

+ 

Starting  date  of  flow  data  for  identification 

only.  Can  be  alphanumeric  data. 

3-10 

or 

QII 

+ 

Flow  data  in  cfs  or  m3/SEC  if  INPUT  (J1.3) 

equals  5.  Remaining  data  starts  in  field  1 

of  succeeding  cards.  (NPER  values.) 

1-10 

QII 

+ 

Flow  data  starts  in  field  1  of  the  second 

card  if  INPUT  (J1.3)  equals  0  (NPER  values). 

*Repeat  IN  cards  for  both  control  points.  Two  sets  of  flow  data  in  turn  for 
each  flood  up  to  the  number  of  floods  prescribed  (NFLOOD  on  J1.7). 
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SUMMARY  OF  INPUT  CAROS 

STREAMFLOW  ROUTING  OPTIMIZATION  PROGRAM  (OPROUT) 


